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ABSTRACT 

We report on multi-object spectroscopy in the red spectral region of 37 candidate star clusters 
in an ~ 8 x 8 arcmin 2 field centered on the giant early-type radio galaxy NGC 1316 (For- 
nax A), the brightest galaxy in the Fornax cluster. Out of this sample, 24 targets are found to 
be genuine star clusters associated with NGC 1316, and 13 targets are Galactic foreground 
stars. For the star cluster sample, we measure a mean heliocentric velocity Vhel = 1698 ± 46 
km s _1 and a velocity dispersion a = 227 ± 33 km s _1 within a galactocentric radius of 24 
kpc. Partly responsible for the velocity dispersion is a significant rotation in the star cluster 
system, with a mean velocity of ~ 175 ± 70 km s _1 along a position angle of ~ 6° ± 18°. 
Using the projected mass estimator and assuming isotropic orbits, the estimated total mass is 
(6.6 ±1.7) x 10 11 M Q within a radius of 24 kpc. The mass is uncertain by about a factor of 
two, depending on the orbital assumptions. The implied Ml Lb ratio is in the range 3-6. Four 
star clusters in our sample are exceptionally luminous (My < —12.3). This means that (1) at 
least this many clusters in NGC 1316 are up to an order of magnitude more luminous than the 
most luminous star cluster in our Galaxy or M 31, and (2) that the S/N ratio of their spectra 
allows us to measure line strengths with good accuracy. By comparing the measured colours 
and equivalent widths of Ha and the Ca II triplet (AA 8498, 8542, 8662 A) absorption lines 
for those bright star clusters in our sample with those of single-burst population models (the 
Bruzual & Chariot 1996 models), we find that they are coeval with an age of 3.0 ± 0.5 Gyr. 
Their metallicities are found to be solar to within ±0.15 dex. We discuss the properties of the 
main body of NGC 1316 and conclude they are consistent with having hosted a major merger 
3 Gyr ago as well. The presence of intermediate-age globular clusters in NGC 1316 shows 
once again that globular clusters with near-solar metallicity do form during galactic merg- 
ers, and, moreover, that they can survive disruption processes taking place during the merger 
(e.g., dynamical friction, tidal disruption), as well as evaporation. In this respect, NGC 1316 
provides a hitherto "missing" evolutionary link between young merger remnants of age ~ 0.5 
Gyr such as NGC 3597, NGC 3921 and NGC 7252 on one side, and older giant ellipticals 
featuring bimodal colour distributions on the other side. 

Keywords: galaxies: individual: NGC 1316 - galaxies: elliptical -galaxies: radio -galaxies: 
interactions - globular clusters: general 



1 INTRODUCTION 1.1 Extragalactic Globular Cluster Systems 

Recent observations with the Hubble Space Telescope (HST) and 
large-field ground-based CCD cameras have caused rapid advances 
t Based on observations obtained at the European Southern Observatory, in our knowledge of the formation and evolution of star cluster sys- 

La Silla, Chile (Observing Programme 60.E-0781). terns of galaxies. Extragalactic star clusters have in recent years es- 

t Electronic mail: goudfroo@stsci.edu tablished themselves as potential tracers of the formation and evo- 
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lution of galaxies. The number of detailed photometric studies is 
increasing steadily, and these studies reveal interesting connections 
between star cluster systems and their host galaxies. 

One important, well-known aspect of star cluster systems 
among galaxies is that the number of star clusters per unit galaxy 
luminosity (named the specific frequency Sn) increases system- 
atically from late-type to early-type galaxies, being ~ 2-3 times 
higher in elliptical (E) than in spiral galaxies of type Sb and later 
(Harris & van den Bergh 1981; Harris 1991). This fact has been 
used as an argument against the scenario for forming ellipticals 
through mergers (e.g., van den Bergh 1995). However, recent obser- 
vations with the HST have led to a wealth of discoveries of young 
star clusters in merging and starburst galaxies. Indeed, it now seems 
that star clusters in general (and globular clusters in particular) 
may form preferentially in high-density regions of starbursts (e.g., 
Meurer et al. 1995), probably from giant molecular clouds whose 
collapse is being triggered by a 100- to 1000-fold increase in gas 
pressure due to supernova and shock heating during starbursts (e.g., 
Elmegreen & Efremov 1997). Galactic gas-rich mergers are known 
to produce the most energetic known starbursts (e.g., Sanders & 
Mirabel 1996), to create galaxy remnants that have surface bright- 
ness profiles consistent with those of "old" giant ellipticals (Wright 
et al. 1990), and to create large numbers of young globular clusters 
(e.g., Whitmore 1999 and references therein). It is possible that the 
higher specific frequency of globular clusters in ellipticals relative 
to that in spirals may be accounted for by secondary populations of 
globular clusters created during gas-rich mergers. 

In this respect, one particularly interesting feature of the star 
cluster systems of many giant ellipticals is the presence of bi- 
modal colour distributions, providing clear evidence for the occur- 
rence of a "second event" in the formation of these systems (e.g., 
Zepf & Ashman 1993; Whitmore et al. 1995; Geisler, Lee & Kim 
1996; Kissler-Patig et al. 1997; Carlson et al. 1998). While such 
a bimodal colour distribution was actually predicted from merger 
models of E galaxy formation by Ashman & Zepf (1992), opin- 



ions about the general nature of the "second event" differ among 
authors (see detailed review by Ashman & Zepf 1998). In par- 
ticular, the "second event" could actually be a series of different 
events that lead to the formation of metal-rich globular clusters. 
It is therefore important to obtain any additional evidence for (or 
against) a merger origin of the star cluster bimodality. One impor- 
tant constraint to any scenario to explain the bimodality is set by 
the colour of the "red" peak in the colour distribution among giant 
ellipticals, which (if interpreted in terms of metallicity) typically 
indicates a near-solar metallicity (e.g., Forbes, Brodie & Grillmair 
1997 and references therein). This metallicity estimate has since 
been confirmed by Keck/LRIS spectroscopy of globular clusters 
in M87 (Cohen, Blakeslee & Rhyzov 1998) and in NGC 1399 
(Kissler-Patig et al. 1998). In order to render the scenario in which 
the second-generation "red" star clusters are formed during galaxy 
mergers generally applicable, it is therefore of great interest to find 
out just what the metallicities of the young, luminous star clusters 
in merger remnants actually are, and whether or not they can be 
found in merger remnants at different times after the last (major) 
merger. 

The current evidence that young star clusters in mergers have 
approximately solar metallicities is still sparse. Photometry of these 
clusters typically does not provide enough information, due to the 
age-metallicity degeneracy of (especially optical) colours and the 
fact that most merger remnants contain extensive dust features with 
appreciable extinction. Progress has recently been made by us- 
ing combinations of near-infrared and optical colours (Maraston 
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Figure 1. Grey-scale plot of B-band CCD image of the central 1'.5 X 1'.5 
of NGC 1316 (taken from Paper II). The plate scale is indicated by the bar 
at the lower left of the image. North is up and East is to the left. 



et al. 2000), but this method is still in its infancy. The only de- 
pendable evidence stems from relatively strong metallic features 
observed in spectra of one bright cluster in the peculiar elliptical 
galaxy NGC 1275 (Zepf et al. 1995a; Brodie et al. 1998) and of 
two such clusters in NGC 7252 (Schweizer & Seitzer 1998), a pro- 
totypical, 0.5 - 1 Gyr old remnant of two merged spiral galaxies 
(Hibbard & Mihos 1995; Schweizer 1998). It is clear that the need 
for more spectroscopic metallicity estimates for young star clusters 
in (other) mergers is imminent. 

Spectroscopy of star clusters around ellipticals can also pro- 
vide one with the kinematics of the cluster system, another es- 
sential source of information for discriminating between differ- 
ent galaxy formation scenarios. For example, spectroscopic stud- 
ies of M87 and NGC 1399 revealed a significantly higher veloc- 
ity dispersion for the globular clusters than for the stars (Mould et 
al. 1990; Brodie and Huchra 1991; Grillmair et al. 1994; Cohen & 
Ryzhov 1997; Minniti et al. 1998; Kissler-Patig & Gebhardt 1998; 
Kissler-Patig et al. 1999). As both M 87 and NGC 1399 are central 
cluster ellipticals, these authors suggested that the globular clus- 
ters were reacting to the gravitational potential of the galaxy clus- 
ter as a whole rathe r than that of the galaxy itself. In NGC 5128 
(Cen A), Hui et al. (1995 ) reported rotation in the globular cluster 
system, though only for the metal-rich clusters which seemed to ro- 
tate along with the well-known dust lane. However, the situation in 
NGC 4472 (Sharpies et al. 1998) and M87 (Kissler-Patig & Geb- 
hardt 1998) seems to be different in that the metal-poor globular 
clusters seem to dominate the rotation. 

The present paper describes new spectroscopic observations 
aimed at deriving kinematics, ages and metallicities for a suitable 
sample of star clusters in the merger remnant NGC 1316. 



1.2 NGC 1316 

The giant early-type galaxy NGC 1316 (= Fornax A = PKS 
0320 - 37 = ESO 357 - G022 = ARP 154) is one of the brightest 
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Table I. Global properties of NGC 1316. 
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SO 1 pec 




r e n 


109" 


CCD94 


CO 


221 km s" 1 


K00 




9.42 


RC3 


foreground 


0.00 


BH84 


(B - V) e « 


0.93 


RC3 


(U - B) cB 


0.47 


RC3 


Vhel 


1760 km s" 1 


NED 


Distance 


22.9 Mpc 


** 


M B T 


-22.38 


* * * 



Notes to Table |l| 

Parameters: r e ff = Effective radius; cro = central velocity dispersion; 
Bj< = Total B magnitude; Ag^ foreground = Absorption in B band due 
to ISM in our Galaxy; (B — V) e g, (U — B) a g = Mean colours inside the 
effective radius; u ne ] = Heliocentric velocity; Mg = Total absolute B 
magnitude, corrected for Galactic absorption. 

References: RC3 = de Vaucouleurs et al. (1991); RSA = Sandage & Tam- 
mann (1987); CCD94 = Caon et al. (1994); K00 = Kuntschner (2000); 
BH84 = Bu rstein & Heiles f!984): NED = NASA/IPAC Extragalac- 
tic Database (http://nedwww.ipac.caltech.edr); * = derived from archival 



HST/WFPC2 images; * * = derived from supernova light curves, see text; 
* * * = derived from B T and A Sj foreground- 

and closest radio galaxies in the sky. It is located in the outskirts of 
the Fornax cluster, at a projected distance of 3? 7 from NGC 1399, 
the central giant elliptical galaxy. The extensive optical observa- 
tions of Schweizer (1980) showed that NGC 1316 is a typical Mor- 
gan D-type galaxy with an elliptical spheroid embedded in an ex- 
tensive, box-shaped envelope. Several features of NGC 1316 estab- 
lish firmly that it is a merger remnant. The outer envelope includes 
several non-concentric arcs, tails and loops that are most likely rem- 
nants of tidal perturbations, while the inner part of the spheroid 
is characterized by a surprisingly high central surface brightness 
and small core radius (and effective radius) for its galaxy luminos- 
ity (Schweizer 1981; Caon, Capaccioli & D'Onofrio 1994). These 
characteristics, together with a velocity dispersion which is signifi- 
cantly lower than that of other elliptical galaxies of similar lumi- 
nosity, cause NGC 1316 to lie far off the fundamental plane of 
early-type galaxies or the Faber-Jackson relation (e.g., D'Onofrio 
et al. 1997). Furthermore, the inner spheroid is marked by non- 
concentric ripples of 0.1-0.2 m ag a mplitude above a best-fitting 
r 1//4 law (see discussion in Sect. 5.2 3). Figure [l] shows a 73-band 
image of the inner 7'5 x 7'5 of NGC 1316 (taken from Goud- 
frooij, Alonso & Minniti 2000) in which some of the features men- 
tioned above can be found. All of these features are consistent with 
NGC 1316 having undergone a relatively recent merger after which 
dynamical relaxation has not yet had time to complete fully. 

As to the adopted distance of NGC 1316, we take advantage 
of the fact that two well-observed type la supernovae (SNela) have 
occurred in NGC 1316 (SN1980N and SN1981D). Using the pre- 
cise distance indicator for SNe la that utilizes the tight relation be- 
tween their light curve shape, luminosity, and colour (Riess et al. 
1998; A. G. Riess, private communication), we arrive at a distance 
of 22.9 (± 0.5) Mpc for NGC 1316, equivalent to (m - M) = 
(31.80 ± 0.05). At this distance, 1" corresponds to 111 pc. Note 
that this sets NGC 1316 and the sub-group of galaxies surrounding 
it slightly behind the core of the Fornax cluster (for which (m— M)o 



= 31.54 ± 0.14; Ferrarese et al. 2000). Global galaxy properties of 
NGC 1316 are listed in Table @. 

This paper is built up as follows. Section ^ describes the sam- 
ple selection. The observations are described in section ^|, while 
section^] deals with the data reduction. The various results are pre- 
sented and discussed in sections ^ and ^j, respectively. Finally, sec- 
tion [7] summarizes the main conclusions of this study. 



2 SELECTION OF STAR CLUSTER CANDIDATES 

The sample of star clus ter ca ndidates was prepared from the tar- 
get lists of Shaya et al. (1996) who used the first-generation Plan- 
etary Camera aboard HST (V,I imaging) and Goudfrooij et al. 
(2000, hereafter referred to as Paper II), who used ground-based 
(B, V, I, J, H, K') imaging supplemented with the HST/WFPC2 
imaging first published by Grillmair et al. ( 1999] ). From these two 
lists of point-like sources around NGC 1316, a first cut was made 
by selecting only targets brighter than V = 22 mag, which is the 
approximate limiting magnitude to obtain a signal-to-noise (S/N) 
within an exposure time of about 4-5 hours on a 4m-class telescope 
that is high enough for reliable velocity measurements. Apart from 
the bright stars that were chosen as pointing checks (see next Sec- 
tion), a further cut was then made by selecting only targets with 
V <; 18 (equivalent to My <; —13.8 at the assumed distance of 
NGC 1316), in order not to miss any possible massive, young or 
intermediate-age star clusters. As to the colour selection criterion, 
we made the following consideration. The colour distribution of 
compact objects in NGC 1316 is very broad (0.8 ^ 73-/ ^ 4.5, 
cf. Paper II), peaking at about B — I — 1.9. In view of the fact 
that the central square arcmin of NGC 1316 is full of dust lanes 
and patches featuring large extinction values which may redden a 
number of star clusters, we considered the full observed range in 
B — I for our target selection. A colour-magnitude diagram for the 
observed targets is shown in Fig. ^. The symbols show the different 
classes of targets, as resulting from the spectral data shown in this 
paper. A histogram of the colour distribution of compact objects in 
the field of NGC 1316 (from Paper II) is added to the figure for 
comparison. 



3 OBSERVATIONS 

Spectroscopy of the star cluster candidates in NGC 1316 was car- 
ried out during November 15-17, 1997, with the ESO 3.5-m New 
Technology Telescope (NTT) equipped with the ESO Multi-Mode 
Instrument (EMMI) in Multi-Object Spectroscopy (MOS) mode. 
The observational setup was identical to the one we described in 
detail in Minniti et al. (1998). Here we give a somewhat shorter 
description. 

The detector was a thin, back-illuminated CCD of type Tek- 
tronix TK2048EB Grade 2, having 2048x2048 sensitive pixels (of 
size 24 /im = 0"268 pixel -1 ). We used grism #4 (2.8 A pixel" 1 , 
blazed at 6500 A). The wavelength range is 5500 - 10000 A. The 
MOS slitlets had a width of l'.'34, resulting in a wavelength res- 
olution of 7.5 A. We chose grism #4 in order to include the Call 
triplet (AA 8498, 8542, 8662 A). The Ca II triplet index is an excel- 
lent metallicity indicator, being independent of age for stellar popu- 
lations older than about 1 Gyr (e.g., Bica & Alloin 1987; Diaz, Ter- 
levich & Terlevich 1989; Garcia- Vargas, Molla & Bressan 1998). 

Once the final target selection was made, it was decided to ro- 
tate the instrument by 90° because that made it possible to include 
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4 DATA REDUCTION AND VELOCITY 
MEASUREMENTS 

The data reduction of the multi-slit spectra was performed in a 
standard manner using the IR AF IM RED.SPECRED package, just 
as described in Minniti et al. (1998). The only noteworthy chal- 
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Figure 2. Optical colour-magnitude diagram for the sample targets. The 
open circles turned out to be foregr ound stars after analysis of the spectra 
presented in this paper (cf. Section p.l| ), and the filled circles turned out 
to be genuine star clusters associated with NGC 1316. A histogram of the 
colour distribution of compact objects in the field of NGC 1316 is shown at 
the bottom for comparison (taken from Paper II). 



all selected targets in only 2 MOS masks (covering 5' x 8' each), 
thus making the most efficient use of telescope time. Before punch- 
ing the slits onto the EMMI MOS plates, the distortion in the focal 
plane of the EMMI instrument has to be derived. To this end, short- 
exposure R-band images were kindly taken by the ESO NTT team 
during the night before our observing run, also using EMMI in the 
rotated position. The (x, y) positions of our targets were measured 
from these images to produce the MOS masks. In each mask, a 
few slitlets were punched at positions corresponding to bright fore- 
ground stars in the field of view, which were used to check the tele- 
scope pointing. Furthermore, a few targets were observed through 
both masks to enable one to check for any systematic offsets be- 
tween radial velocities derived using the different masks. One field 
was centered about 100" north of the centre of NGC 1316, and 
observed for a total exposure time of 14,400 s, and the other field 
(centered about 90" south of NGC 1316, and including star cluster 
candidates near the galaxy centre which, due to the higher galaxy 
surface brightness, require a longer exposure time) was exposed 
for a total of 17,900 s. Individual MOS exposures were typically of 
1800 s duration, interleaved with HeAr lamp exposures and short 
imaging exposures to correct for small pointing errors due to flex- 
ure and differential atmospheric refraction. The latter corrections 
were always less than 0'.'3, and the air mass was always kept below 
1.3. The weather was photometric throughout the observing run, 
and the seeing varied from (X'6 to 0V8 as judged from the direct 
images taken during the night. 



lenge was a proper subtraction of the local sky background from 
the targets that are close to the galaxy centre, where the galaxy 
background light is changing rather rapidly with radius. However, 
this problem was foreseen: We punched longer slits for those tar- 
gets near the galaxy centre, enabling us to make reliable third- 
order polynomial fits to the galaxy background level as a func- 
tion of position along the slit, which usually provided good results. 
The only two targets where special care was necessary were tar- 
gets #112 and #114 (= knots D and E in Schweizer 1980, re- 
spectively) (cf. Table g), since they are situated only 14'.'5 and 
33'.'9 from the centre, respectively, within a filament system of ion- 
ized gas (cf. Schweizer 1980; Mackie & Fabbiano 1998). This in- 
troduces the additional problem of the presence of radially vari- 
able emission lines. This problem was dealt with as follows. We 
first subtracted the emission-line-free background light with the 
usual third-order polynomial fit to the background light, but ex- 
plicitly excluding the wavelength regions around the emission lines 
[O I] AA6300, 6363, [N II] AA6548, 6583, [S II] AA6716, 6731, and 
[S III] AA9069, 9531 (redshifted according to the systemic velocity 
of NGC 13 16) in the fitting procedure. The next (and final) step was 
to fit a polynomial to the positions and intensity of those emission 
lines in the background-subtracted spectra as a function of posi- 
tion along the slit, as if they represented a separate system of tel- 
luric lines. The background in between the emission lines associ- 
ated with NGC 1316 was median filtered during the creation of this 
second "emission-line sky" spectrum to avoid any further degrada- 
tion of the S/N ratio of the final object spectrum. This procedure 
proved successful in the case of target #114, since its resulting 
spectrum (cf. Fig. [§]) is free of any residual signs of [N II] A6583 
emission. However, the resulting spectrum of target #112 still con- 
tains residual Ha and [N II] A6583 emission lines, presumably due 
to its location on a compact, localized emission-line filament which 
is not present on the sky along the slitlet. Hence, target #112 was 
excluded from the discussion on line strength indices in Section 



After wavelength calibration, the spectra were corrected for 
atmospheric extinction and flux calibrated usi ng a spectrum of th e 
spectrophotometric standard star LTT 1020 (Hamuy et al. 1994), 
which was taken through one of the MOS slitlets at an airmass sim- 
ilar to those of the NGC 1316 exposures. 

The final spectra were cross correlated with s tellar template 
spectra using the method of Tonry & Davis ( 1979] ), incorporated 
in the FXCOR task within IRAF. As template spectra we used both 
the radial velocity standard HR 1285 (a KOlII star, taken though 
a MOS slitlet), and two bright foreground stars that were already 
observed as pointing reference stars. The spectral regions covered 
by molecular absorption lines due to the Earth atmosphere were 
explicitly excluded in the cross correlations. The zero point of the 
radial velocities was good to 50 km s _1 RMS, as judged from the 
cross correlations of the different templates against one another. 



5 RESULTS 
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Table 2. Star clusters and foreground stars: Positions, velocities, photometry, and classifications. 
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RA 
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(J2000) 
(3) 


fhel 

(kms" 1 ) 
(4) 


Av 
(kms" 1 ) 
(5) 
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(V) 
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3 


36 


50 


21.37 


1.38 


3.38 


Star 


1 13 


03 22 40.73 


—37 09 24.1 


291 


334 


107 


22.17 


1.52 


3.33 


Star 


1 17 


03 22 46.55 


—37 12 27.7 


—21 


46 


64 


20.82 


1.31 


3.54 


Star 


118 


03 22 47.56 


-37 13 18.0 


15 


12 


47 


20.68 


1.17 


2.55 


Star 


1Z4 


m n 

U3 Zz jj.jo 


— 3 1 lz jZ.Z 


OA 
ZU 


z 


JO 


on AQ 
ZU.Oo 


1 A q 
1 .48 


1/11 

3.43 


C + nr 

Mar 


125 


03 22 56.64 


—37 11 36.8 


253 


13 


60 


20.59 


0.54 


1.17 


Star 


126 


03 22 58 2 3 


— 37 1 1 38 ? 


7 


o 


30 


17.63 


1.04 


2.21 


Rpfprpnrp <;tnr 


202 


03 22 25.02 


-37 12 45.1 


-6 


15 


72 


19.27 


1.56 




Star 


206 


03 22 32.11 


-37 13 07.3 


35 


9 


173 


22.02 


1.69 


4.31 


Star 


209 


03 22 37.04 


-37 14 33.6 


142 


1 


29 


18.26 


0.73 


1.56 


Reference star 


213 


03 22 42.14 


-37 11 14.6 


430 


5 


73 


21.27 


0.67 


1.98 


Star 


214 


03 22 44.33 


-37 12 44.9 


2 


4 


49 


16.26 


0.28 


1.98 


Reference star 


218 


03 22 51.17 


-37 14 20.4 


165 


23 


65 


20.94 


0.39 


0.91 


Star 


219 


03 22 52.26 


-37 13 33.9 


126 


19 


172 


22.36 


1.26 


3.49 


Star 


220 


03 22 55.10 


-37 15 00.0 


240 


57 


82 


20.75 


0.46 


1.03 


Star 



Notes to Table Column (1) contains the target ID; columns (2) and (3) contain the J2000.0 equatorial coordinates. Units of RA 
are hours, minutes, and seconds; units of DEC are degrees, arcminutes, and arcseconds. For comparison, the J2000 coordinates of the 
nucleus of NGC 1316 are (RA, DEC) = (03 d 22 m 41?710, -37° 12'2#'41). Column (4) lists the heliocentric radial velocities; columns 
(5) and (6) list the velocity uncertainties (At) gives the formal error from the cross-correlation; e v gives the maximum difference of 
velocities derived using the different templates). Columns (7), (8) and (9) list the magnitudes and colours, respectively, and column 
(10) lists the object classification assigned by us. "Reference" stars designate stars that were used for pointing checks, cf. Section Bl 

5.1 Star cluster system membership 

Table ^ lists the observed target list of star cluster candidates in 
NGC 1316. For each target we give an ID number, the J2000 equa- 
torial coordinates, the heliocentric radial velocity and its uncertain- 
ties, the photometry in V, B — V and B — I (taken from Paper II), 
and the assigned nature of the target (star/ globular). The table has 
been sorted to separate stars and globular clusters for convenience. 

A total of 17 targets are found to be Galactic foreground stars 
© 0000 RAS, MNRAS 000, 000-000 



(see Table |2j). Four of these were deliberately chosen as pointing 
alignment check (cf. Section ||); the other 13 stars turned out to be 
failed star cluster candidates. The remaining 24 targets are genuine 
star clusters associated with NGC 1316. This represents an overall 
"success rate" of 65%, although this is strongly a function of colour 
(e.g., it is 80% in the range 1.2 < B — I < 2.4 in which old, non- 
reddened clusters are expected). Figure |3| depicts the dependence 
of the heliocentric radial velocities on the V magnitudes and B — I 
colours for the foreground stars and the NGC 1316 star clusters. 
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Figure 3. Dependence of heliocentric radial velocity on V magnitude (left) and B — I colour (right) for the foreground stars (open circles) and star clusters 
(filled circles) of the sample. 
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Figure 4. Heliocentric radial velocity distribution of the 24 star clusters in 
NGC 1316. The dashed line represents the mean velocity of the star clusters, 
and the dotted line represents the systemic velocity of the integrated light 
of NGC 1316. 



Note the clear separation in radial velocity between the foreground 
stars and the star clusters, and the rather narrow distribution of B—I 
colours of the genuine star clusters. Only two of the selected "red" 
targets turned out to be star clusters. All other targets with B — I> 
2.1 are foreground M stars. 



5.2 Analysis of the star cluster velocity field 



5.2.1 Velocity distribution of the star clusters 

Figure ^ shows a histogram of the heliocentric velocities of the 
star clusters in our sample. The velocity distribution exhibits two 
peaks: one on either side of the systematic velocity of the galaxy. 
Fitting a Gaussian to this distribution yields a mean heliocentric 
velocity Vhel = 1698 ± 48 km s _1 , similar to the systemic velocity 
of NGC 1316 (cf. Table |l]). The formal velocity dispersion of the 
full sample is a — 227 ± 33 km s _1 . However, the double-peaked 
velocity distribution signals caution to be taken before assigning 
a physical meaning to this value of a. Our interpretation of the 
double-peaked velocity distribution is a rotation of the star cluster 
system, cf. below. 



5.2.2 Spatial distribution of the velocities 

In Figure ^ we depict the relation of the star cluster velocities with 
their position angle (hereafter PA; measured from the nucleus of 
NGC 1316, North through East), as a function of projected galac- 
tocentric radius. Apart from a few outliers (e.g., the star cluster at 
PA ~ 305°, Mhel ~ 1250 km s _1 ), the distribution of the velocities 
in this plot has a periodic appearance. Interpreting this as rotation 
of the star cluster sample, we fitted a fiat rotation curve to the he- 
liocentric velocities using the parametric function 

Vhel = Vsys + vo cos(6> — f?o ) 

where v sys is the systemic velocity, vo is the rotation velocity, 9 is 
the PA of the star cluster, and 8q is the PA of the kinematic ma- 
jor axis of the star cluster system. The three free parameters were 
solved for numerically. The uncertainty of the fit is provided by 
the standard covariance matrix. The uncertainties of the individ- 
ual velocities were also accounted for. However, the uncertainty 
of the derived rotation velocity is dominated by the high velocity 
dispersion of the star cluster system. The latter was derived after 
subtracting the best-fitting flat rotation curve. The fitting procedure 
was done twice: once by leaving v sys a free parameter, and once by 
fixing v sys to be the systemic velocity of the stellar component of 
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Figure 5. Radial velocities of the star clusters vs. position angle measured 
from the nucleus of NGC 1316 (North through East). The symbol size scales 
linearly with galactocentric radius. The dashed line represents the mean 
radial velocity of the star clusters, the dotted line represents the systematic 
velocity of NGC 1316, and the solid line represen ts a maximum-likelihood 
cosine fit to the star cluster velocities (see Section |5.2] for details). 



Table 3. Results of kinematic fitting of star cluster velocities. 



Rotation 


Velocity 


Position 


Mean 


velocity 


dispersion 


angle 


velocity 


(kms" 1 ) 


(kms" 1 ) 


(degrees) 


(kms" 1 ) 


(1) 


(2) 


(3) 


(4) 


179 ± 68 


202 ± 33 


1 ± 17 


1698 ±48 


172 ± 71 


212 ±33 


10 ± 19 


1760* 



NOTE — * Fixed mean velocity to be NGC 1316's systematic velocity. 

NGC 1316, 1760 km s~ . The results of this parametric fitting of 
the velocity field are listed in Table ^. 

No obvious trends of the rotation velocity with galactocentric 
radius are apparent from Figure ||. This is shown more directly in 
Figure ^, where we plot radial velocity vs. galactocentric radius for 
the star clusters. One possibly relevant feature in this plot is that the 
two innermost star clusters (within 30" from the centre) both have 
a radial velocity that is ~450 km s _1 lower than the systematic 
velocity of NGC 1316. This difference is (in an absolute sense) 
~ 4 a higher than the rotation velocity derived from the full star 
cluster sample. This might indicate either a significant amount of 
radial motion for those innermost star clusters, or the presence of 
differential rotation in the inner regions. 



5.2.3 Comparison with the kinematics of stars and planetary 
nebulae 



Arnaboldi et al. ( 1998 ) studied the dynamics in the outer regions of 
NGC 1316 using MOS spectra of planetary nebulae in NGC 1316. 
They find that the planetary nebulae in NGC 1316 rotate along 
PAp N = 80°±30° with a rotation velocity v PN = (140 ± 60) 
km s _1 at a galactocentric radius of 4', and a velocity dispersion 
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Figure 6. Radial distribution of the velocities of star clusters in NGC 1316. 
Filled symbols represent star clusters on the north side of the apparent major 
axis of the stellar body NGC 1316; open symbols represent the star clusters 
on the south side of it. Dashed and dotted lines as in Figure |^. 



crpN = (175 ± 18) km s _1 . Note that while the rotation velocity 
and the velocity dispersion of their planetary nebula sample agree 
with those of our star cluster sample to within the combined 1 a 
error, the PA of the kinematic major axis of the planetary nebula 
sample is different from that of the star cluster sample by ~ 70° 
(although this is formally only a 2.5 a result). 

As to the kinematics of the integrated stellar light of 
NGC 1316, the PA of the kinematic major axis appears to be con- 
sistent with that of the apparent major axis of the isophotes in the 
inner ~60", PA» ~ 58°(Arnaboldi et al. 1998; Bosma, Smith & 
Wellington 1985). Both Arnaboldi et al. and Bosma et al. report 
that the stellar velocity field in the inner ~ 30" is consistent with 
NGC 1316 being an oblate, axisymmetric spheroid flattened by ro- 
tation and having isotropic velocity dispersion. However, irregu- 
lar kinematic features are found at larger distances along both the 
major and minor axes: A jump of ~ 30 km s _1 in radial velocity 
along the major axis is found at the position of a ripple (Bosma 
et al. 1985), and similar peculiar velocities and velocity dispersion 
values are also found along the minor axis (Arnaboldi et al. 1998). 
These features occur at positions where prominent shells or rip- 
ples are present (as checked on our direct CCD images). We argue 
that these features do not affect the global dynamics significantly, 
but are associated with the recent merger history of NGC 1316. In- 
deed, several authors have argued in the past that shells and ripples 
such as those found in NGC 1316 can be formed during a galaxy 
merger (e.g., Quinn 1984; Hernquist & Quinn 1988; Hernquist & 
Spergel 1992) or during a weak interaction (Thompson 1991). The 
shells and ripples reflect the various stages of the wrapping of the 
stars from the smaller galaxy around the centre of the potential. 
Peculiar velocities are indeed expected to be associated with the 
ripples, and the observed velocity jumps of 30 - 50 km s _1 are in 
agreement with the theoretical expectations (cf. Quinn 1984). 
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5.2.4 Mass estimates 



5.3 Star cluster ages and metallicities 



Using the globular clusters in our sample as test particles within 
the underlying galaxy potential well, one can estimate a dynamical 
mass for the galaxy and derive its mass-to-light ratio for the inner 
217" (= 24 kpc, the projected distance of our outermost cluster). 
We estimate the mass using all 24 velocities from globular clusters 
and applying the virial mass estimator and projected mass estima- 
tor (see, e,g., Bahcall & Tremaine 1981; Heisler et al. 1985). No 
assumption is made concerning the systemic velocity which was 
allowed to vary over a wide range, and masses are computed for 
various assumptions concerning the cluster orbits (radial, tangen- 
tial, mixed, or isotropic), following Puzia et al. (2000). The results 
are displayed in Figure M. 

Note that the minimum masses are derived for systemic ve- 
locities close to 1760 km s - , the heliocentric velocity of the in- 
tegrated light of the galaxy. The minimum masses vary between 
4.4 x 10 11 M Q and 1.3 x 10 12 M depending on the orbital 
assumptions. Each measurement has an internal uncertainty of 
~ 30% due to the limited spatial sampling, velocity errors, and 
sample size, and could rise by ~ 20% by assuming a systematic 
velocity that is off by 100 km s _1 . The masses are also dependent 
on the assumed distance. For instance, if NGC 1316 would be at 
the distance of the central part of the Fornax galaxy cluster as de- 
termined by a compilation of Cepheids, Surface Brightness fluctu- 
ations and Planetary Nebulae luminosity functions in other Fornax 
cluster galaxies ((m-M)o = 31.30±0.04, Kavelaars et al. 2000), 
it would reduce the masses by ~ 30%. The dynamical masses de- 
termined using the star clusters as test particles are consistent with 
those determined by Amaboldi et al. (1998) using planetary nebu- 
lae, taking the difference in the assumed distances into account. 



5.2.5 Mass-to-light ratios 

As a working hypothesis to derive mass-to-light (MIL) ratios for 
the stellar component of NGC 1316, we will use the mass derived 
for isotropic orbits at the measured systemic velocity: (6.6± 1.7) x 
10 11 M . 

Total galaxy luminosities in B, I, and H were measured in 
a circular aperture of 220" radius (the galactocentric radius of the 
outermost star cluster) using our optical and near-infrared surface 
photometry of the galaxy (presented in Paper II). 

The results are listed in Table The table also lists total lu- 
minosities for the integrated light of NGC 1316 after dimming to 
an age of 14 Gyr. To estimate the latter, we assumed that the stars 
within a 20" radius from the centre are 2 Gyr old (which is the 
luminosity-averaged age of the central 2 x 4", spectroscopically 
measured by Kuntschner 2000). Outside the 20" radius, we as- 
sume that the bulk of the stars are already old, since the colours 
and colour gradients in that radial range are consistent with those 
of "normal" giant elliptical galaxies (Goudfrooij, in preparation; 
Schweizer 1980). The actual dimming factors of the composite 
population were calculated using the population synthesis models 
of Bruzual & Chariot (1996). The resulting age dimming is only 
of order 0.1 mag in each band, both for Salpeter (1955) and Scalo 
(1986) stellar initial mass functions (IMFs). While this calculation 
is only a rough estimate, it should not be off by any significant 
amount unless the stellar population of the whole galaxy would — 
unexpectedly — be only a few Gyr old. However, we note that the 
resulting MIL ratios are rather low (e.g., MILb ~ 4.5), which ac- 
tually does indicate that the mean stellar population in NGC 1316 
may be younger than usual for giant elliptical galaxies. 



During the early 1990's, a number of different groups have devel- 
oped powerful evolutionary population synthesis models employ- 
ing different stellar metallicities (Worthey 1994; Vazdekis et al. 
1996; Bruzual & Chariot 1996, hereafter BC96; Bressan, Chiosi 
& Tantalo 1996; Maraston 1998). Hence, the determination of ages 
and metallicities of single-burst populations such as star clusters is 
now relatively straightforward in principle. However, the determi- 
nation of ages and metallicities of populations of non-solar metal- 
licity is hampered by the lack of model spectra of sufficient spectral 
resolution for metallicities other than solar. The model spectra of 
non-solar metallicity that are used in most spectrophotometric pop- 
ulation synthesis models are based on the Kurucz (1995, private 
communication) atmospheres, which have an intrinsic wavelength 
resolution of 20 A in the optical region. On the other hand, this is 
not a significant concern for the spectral indices we chose for de- 
riving ages and metallicities (Ha and Ca II, see below), since they 
are well resolved in the model spectra. 



5.3.1 Line strength measurements 

To determine star cluster ages and metallicities, the observed spec- 
tra were convolved with a Gaussian kernel in order to match the 
wavelength resolution of the model spectra. Line strengths were 
then measured from these smoothed observed spectra, and com- 
pared with those measured in an identical manner from a set of 
model spectra. To try to achieve a suitable break of the age- 
metallicity degeneracy using lines in our wavelength coverage, we 
measured Ha (being primarily age sensitive) and the Ca II triplet, 
the equivalent width (hereafter EW) of which is controlled by the 
stellar metallicity for stellar populations older than 1 Gyr (Bica 
& Alloin 1987; Diaz et al. 1989; Garcia- Vargas et al. 1998). The 
line str engths were measured as described in Goudfrooij & Em- 
selle m (|1996|) . The errors were quantitatively estimated following 
Rich ( |1988| ). Table ^ lists the measured equivalent widths for the 
three bright NGC 1316 star clusters with ID's 103, 114, and 210. 
The spectra of the fainter clusters had a S/N too low for reliable 
line strength measurements. Figure || shows the flux-calibrated, ex- 
tracted spectra of these star clusters. Their spatial position within 
NGC 1316 is illustrated in Figure ^ which depicts the l/-band 
"residual image" constructed from a V-band image after subtrac- 
tion of an elliptical isophotal model image (details on this are fur- 
ther discussed in Paper II). 

The adopted passband for EW measurements of Ha was 24 A 
wide. As to the Call triplet measurements, we adopted the CaT 
index 

EW(CaT) ee EW(CaII 8542) + EW(CaII 8662) 



as defined by Diaz et al. (1989). Their definition is well suited for 



the resolution of the model spectra. 



5.3.2 Age and metallicity estimates 

Figures 0and[ll| illustrate the method used to derive cluster ages 
and metallicities from the measured equivalent widths. Figure [Ul 
shows a diagram of EW(Ha) vs. EW(CaT), in which the drawn 
lines represent the age-metallicity grid as measured from the model 
spectra of different ages and metallicities ( BC96) . We used the 
model spectra of Lejeune, Cuisinier & Buser ( 1991) for our analy- 
sis (see BC96 for details). Both Salpeter (1955) and Scalo (1986) 
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Figure 7. Mass estimates for NGC 1316 out to a radius of r < 217" or R < 24 kpc, as a function of systematic velocity. Different orbit characteristics 
of the globular clusters were assumed for the mass estimate using the Projected Mass Estimator (PME): PME; so assumes isotropic orbits, PME rad assumes 
radial orbits, PME m i x assumes mixed orbits, and PMEtan assumes tangential orbits (for details see Bahcall & Tremaine 1981). VME shows the results for 
the Virial Mass Estimator. The dashed dotted lines show the mean and median velocity of our globular cluster sample; the systemic velocity of NGC 1316 is 
1760 km s" 1 . 

[t] 

Table 4. Mass-to-light ratios for NGC 1316 within 24 kpc radius. 



Filter 


Luminosity (Lq) 


MIL ratio 


Corrected luminosity (Lq) 


corrected MIL ratio 


B 


1.58 X 10 11 


4.2 


1.43 x 10 11 


4.6 


1 


2.70 x 10 11 


2.4 


2.46 X 10 11 


2.7 


H 


7.87 x 10 11 


0.8 


7.18 x 10 11 


0.9 



Notes to TaWe^. 

Luminosities are derived assuming Mg^ = 5.48, Mj,q = 4.07, and Mjj,q = 3.39 (Allen 1973; 
Campins, Rieke & Lebofsky 1985). The listed MIL ratios have been calculated using isotropic 
orbits; it can vary by a factor of 2 depending on the orbit assumptions made 
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Figure 9. Negative of y-band residual image of NGC 1316 (taken from Paper II), centered 47" west of the nucleus. White areas (i.e., negative intensity) 
indicate galaxy light being absorbed by dust. The black elliptical structure (plus white "ring") in the central regio n of NGC 1316 is a CCD defect caused by 
saturation. North is up and east is to the left. The three bright star clusters with ID's 103, 1 14, and 210 (cf. Section |5.3|) are indicated, as is the plate scale. 



Table 5. Equivalent widths, ages, metallicities, and absolute magnitudes of star clusters in NGC 1316 





Ha 


Can 8542 


Can 8662 


CaT 


Log Age a 




M v h 


Cluster ID 


(A) 


(A) 


(A) 


(A) 


(yr) 


(dex) 


(mag) 


103 


1.44 ±0.08 


5.17 ±0.20 


3.55 ±0.16 


8.72 ± 0.26 


9.45 ±0.13 


+0.03 ± 0.04 


-12.5 ± 0.1 


114 


1.47 ±0.06 


4.64 ±0.12 


3.57 ±0.10 


8.21 ±0.16 


9.46 ±0.11 


-0.03 ± 0.03 


-13.0 ±0.1 


210 


1.50 ±0.10 


4.62 ± 0.22 


3.09 ±0.18 


7.71 ±0.28 


9.48 ±0.17 


-0.13 ±0.06 


-12.1 ±0.1 



a Ages and [Z/Zq] values estimated from location of clusters in the Ha vs. CaT diagram depicted in Fig. 
b Derived from V magnitudes in Tablet! assuming m — M = 31.8 ± 0.1. 



stellar initial mass functions (IMFs) were considered. Since the re- 
sults were indistinguishable for our purposes, we only plot the re- 
sults using the Salpeter IMF in the diagrams. Note the diagnos- 
tic power of the EW(Hq) vs. EW(CaT) diagram in breaking the 
age-metallicity degeneracy, especially for stellar populations with 
a (luminosity-weighted) age of 1-3 Gyr. This diagram yields clus- 
ter ages and metallicities encompassing remarkably narrow ranges: 
2.8 & Age <> 3.0 Gyr, and -0.13 <> [Z/Z Q ] £ 0.03 (cf. Table |). 

These age- and metallicity estimates are compared to clus- 
ter colours in Figure [llj The two top panels show the evolution 
of EW(Hq) and EW(CaT) as a function of metallicity, as before. 
The two other panels in Figure |TT] show the evolution of B — V 
and B — I in the Kron-Cousins system, again using the models 
of BC96. For comparison with the models, Figure |ll] also shows 
the measured EWs and colours — marked by horizontal lines — of 
the three bright NGC 1316 star clusters. Note that for about solar 
metallicity (as derived from the CaT measurements), the star clus- 
ter colours indicate Log (Age) ~ 9.45 (i.e., Age ~ 3 Gyr), which is 



completely consistent with those derived from EW(Ha) to within 
the uncertainties. 

Note that the Ha and Ca T equivalent widths of the star clus- 
ters by themselves also allow for an alternative cluster age of ~ 10 7 
yr (cf. Figure |ll|). However, comparison of the colours of the clus- 
ters with those of the BC96 models shows that this alternative 
solution would require those clusters to be highly reddened (by 
0.7 <, E(B-V) <, 0.8). Since two of these three clusters (#103 
and #210) are located well away from the dust features near the 
galaxy centre (cf . Figure ^ and since the spectral continuum shape 
of the three clusters are extremely similar to one another (cf. Fig- 
ure this possibility is extremely unlikely (in fact, the similarity 
between the spectra and colours of the three clusters argues that 
cluster #114 is located in front of the dust patches onto which it is 
superposed on the sky). Hence, we assume in the following that the 
bright star clusters are ~ 3 Gyr old. Obtaining additional high-S/N 
blue-to- visual spectra of these clusters (including the higher-order 
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Figure 8. Flux-calibrated spectra of three bright star clusters in NGC 1316 
(ID's 103, 114, and 210). Atmospheric 2 features (6858 - 6934 A, 7580 - 
7690 A) have been interpolated over. To diminish the influence of noise, the 
spectra have been smoothed with a boxcar kernel of 1 or 2 pixels, depending 
on the S/N ratio of their continuum. Main absorption features are identified 
in the top panel. Note the relatively strong metal features (compare with 
e.g., NGC 1783 and NGC 1978 in Fig. 6b of Bica & Alloin 1987). 



B aimer lines) would certainly enable one to discriminate conclu- 
sively between the two ages described above. 



6 DISCUSSION 

In this section we discuss various issues regarding the physical na- 
ture of the observed star clusters in NGC 1316, and the relevance 
of the age and metallicity of these clusters to ( i ) observed proper- 
ties of the main body of NGC 1316, and (ii) the nature of bimodal 
colour distributions of globular cluster systems in "normal" giant 
elliptical galaxies. 



6.1 Nature of intermediate-age clusters in NGC 1316 

6.1.1 Cluster luminosities and masses 

Given the distance of NGC 1316 ((m-M)o = 31.8), the luminosi- 
ties of the three bright star clusters are high (My = —13.0 for 
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Figure 10. EW(Ha) vs. EW(CaT) for clusters 103, 1 14, and 210 (filled cir- 
cles). A grid of isochrones (solid lines) and iso-metallicity contours (dotted 
lines) based on the Bruzual & Chariot (1996) models is superposed. From 
this diagram, the cluster metallicities appear to be close to solar and the ages 
~ 3 Gyr. 



ID#114, -12.5 for ID#103, and -12.1 for ID#210, cf. Table g). 
For comparison, the most luminous (by far) globular cluster in our 
Galaxy, wCen, has M v = -10.2 (e.g., Meylan et al. 1995), and 
Gl, the most luminous cluster in M 31, has My = —10.85 (Reedet 
al. 1994). According to the BC96 models, a 3 Gyr old cluster with a 
Salpeter IMF fades by only AMy — 1.5 mag over the next 1 1 Gyr 
(or by 1 .4 mag with a Scalo IMF). Hence, at an age of 14 Gyr, these 
three clusters should have absolute magnitudes (My ) of (at least) 
— 11.5, —11.0, and —10.6, i.e., ~2.8, 1.7, and 1.2 times as lumi- 
nous as u> Cen, respe ctively. Assuming Ml Lv = 4.1 as for uj Cen 
(Vleylan et al. 1995), this indicates impressive cluster masses of 



1.4 x 10' M© for #114, 8.5 x 10 6 M for #103, and 6.0 x 10 6 M Q 
for #210, after evolving to an "old" age of 14 Gyr. These masses 
are similar to that of Gl in M 31, the most massive globular clus- 
ter known from kinematic measurements ((1.4 ± 0.4) x 10 7 Mq, 
Meylan et al. 1997, 2000). 



6. 1.2 Open or globular clusters J 



Grillmair et al. (1999) studied the star cluster system of NGC 1316 
using HST/WFPC2 images, whereby they excluded the data within 
a ~ 70" galactocentric radius to avoid the dusty areas. They found 
a luminosity function (LF) of the star clusters that does not resem- 
ble the log-normal globular cluster LF observed in almost all glob- 
ular cluster systems of "normal" (presumed old) elliptical galax- 
ies (see, e.g., Ashman & Zepf 1998). Instead, the LF seemed to 
increase all the way to the faintness limits of their photometry. 
This LF behaviour has been seen in other merger remnants as well, 
NGC 4038/4039 ( Whitmore & Schwe izer 1995; Whitmore et al. 
199 9), NGC 3597 dCarlson et al. 1999b and NGC 3256 (Z epf et al. 
1999) being three well-studied cases. In view of the fact that the 
open cluster population in our Galaxy shows an LF whose shape 
is similar to that of the star clusters in NGC 1316, Grillmair et al. 



(1999) argue that these objects may be young open clusters rather 
t han g lobular clusters (an argument used as well by van den Bergh 
(|l995|) for the case of NGC 4038/4039). However, the new evi- 



dence provided by the much older age of the bright star clusters 
in NGC 1316 (with respect to those in the young merger remnants 
mentioned above) in combination with the high cluster luminosities 
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Figure 11. Time evolution of spectral line indices and colour indices of 
single-burst stellar populations (Brazual & Chariot 1996). Model curves 
are plotted for a Salpeter (1955) IMF and the following metallicities: 0.2 
solar (dot-short-dashed lines), 0.4 solar (short-long-dashed lines), 1.0 solar 
(solid lines), and 2.5 solar (dot-long-dashed lines). The models are com- 
pared with observations (horizontal lines) of star clusters 103 (short-dashed 
line), 114 (dotted line) and 210 (long-dashed line). Equivalent widths of Ha 
and Ca T are in units of A, while B— V and B— I are in units of magnitude. 
The observed line strengths and the colour indices of the star clusters are 
consistent with them having ~ solar metallicity and an age of ~ 3 Gyr. 



and masses renders the nature of these objects to be genuine glob- 
ular clusters beyond any reasonable doubt, at least in NGC 1316. 

As to the apparen t disag reement between our conclusion and 
that of Grillmair et al. ( 19991), we wish to note that only one of the 



three clusters for which we have spectroscopic age measurements 
was covered by the WFPC2 imaging of Grillmair et al., and that the 
size distribution of clusters brighter than / ~ 23.5 is indistinguish- 
able from that in NGC 1399, an old elliptical galaxy (cf. Figure 12 
of Grillmair et al.). On the other hand, we consider it quite likely 
that the faint end of the observed cluster system in NGC 1316 may 
be partly populated by open clusters, in agreement with Grillmair 
et al.'s conclusion. In fact, these are quite likely candidates of star 
clusters that will be destroyed over the next 10-12 Gyr as a result 
of evaporation and/or tidal shocking (see, e.g., Meylan & Heggie 
1997). Gradual removal of these objects from the cluster system 
would conceivably leave behind a log-normal LF as observed in 
cluster systems of "normal" giant ellipticals. This issue will be dis- 
cussed further in Paper II. 



6.1.3 Fate of young and intermediate-age clusters 

Extraordinarily luminous, massive, and metal-rich globular clus- 
ters such as the ones found here are now known to be quite com- 
mon in young merger remnant galaxies, e.g., NGC 1275 (Zepf et 
al. 1995a, Brodie et al. 1998), NGC 7252 (Schweizer & Seitzer 
1998), and NGC 3597 (Carlson et al. 1999), all of which have clus- 
ter ages of order 300-500 Myr. The fate of these massive clusters 
has been an issue of some debate. From spectroscopy of luminous 
clusters in NGC 1275, Brodie et al. (1998) found H 7 and US EWs 
that are somewhat larger than those of any BC96 models that use 
Salpeter or Scalo IMFs. They argued that the large Balmer EWs 
can only be brought into agreement with the BC96 models by as- 
suming a flatter IMF, which they mimic by introducing a low-mass 
cutoff of ~ 2 M Q and a high-mass cutoff of ~ 3 M to the IMF. 
With such a small or absent low-mass stellar component, they ar- 
gue that the young clusters should fade away in only ~ 1 10 9 yr. 
On the other hand, the ~ 500 Myr young clusters in NGC 7252 
show Balmer line strengths that are well reproduced by the BC96 
models that use the standard Salpeter or Scalo IMF (Schweizer & 
Seitzer 1998). Again, the presence of similarly massive, but ~2.5 
Gyr older, globular clusters in NGC 1316 shows that such lumi- 
nous clusters can survive the dynamically hot period during the 
first few Gyr after a galaxy merger. We suggest that these high 
masses be checked through velocity dispers ion m easurements such 
as those performed by Ho & Filippenko (1990 for the brightest 



young cluster in NGC 1569. One should be able to perform such 
measurements with 8- 10 m class telescopes, at least for clusters 
#114 and #103. If the high masses are confirmed for these 3 Gyr 
old clusters in NGC 1316, it would provide strong evidence in sup- 
port of the view that mergers of (gas-rich) galaxies can form very 
massive clusters, with masses similar to or larger than the largest 
ones present in our Galaxy or M 3 1 . 



6.2 A formation scenario for NGC 1316 

How does the ~ 3 Gyr age of the bright clusters in NGC 1316 re- 
late to the properties of its main (stellar and gaseous) body? While 
a multitude of features of NGC 1316 firmly establish its being a 
merger remnant (e.g., tidal tails, shells and ripples, small effective 
radius; see Introduction), the complicated structure of NGC 1316 
has hampered an accurate reconstruction of its merger history. 

However, recent numerical simulations of major disk-disk 
mergers (e.g., Hibbard & Mihos 1995; Bames & Hernquist 1996) 
reveal several features that are very similar to those in NGC 1316. 
In particular, we would like to point out that several properties of 
NGC 1316 are remarkably consistent with it being an evolved rem- 
nant of a merger event similar to the one which formed NGC 7252 
(0.5- 1 Gyr ago, cf. Schweizer & Seitzer 1998). These consisten- 
cies go far beyond the similarity of the masses and metallicities of 
their respective bright globular clusters. 

To set the stage for this comparison, we mention the results 
of Hibbard et al. (1994)'s optical, X-ray, and Hi observations of 
NGC 7252. H I exists in large amounts in its tidal tails. However, 
H I is not detected in its main stellar body. X-ray emission from hot 
(10 6 -10 7 K) gas is detected from the main stellar body, with a lu- 
minosity that is higher than the typical X-ray luminosity of a spiral 
galaxy of the same optical luminosity as NGC 7252. To explain the 
H I observations, Hibbard & Mihos (1995) performed iV-body sim- 
ulations of NGC 7252. They successfully modeled the H I velocity 
reversal in the tails as gaseous material from the pre-merger disks 
that has reached its orbital turnaround point and is falling back to- 
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wards the central regions of the merger remnant. Furthermore, their 
best-fitting numerical model for NGC 7252 raised several issues 
relevant to the current comparison with NGC 1316, which are dis- 
cussed below. 

(i) The fate of cold gas in the progenitor disks. Extrapolating 
from the current H I content of NGC 7252's tidal tails, Hibbard & 
Mihos (1995) found that several 10 9 Mq of atomic gas had already 
flown inward since the merger took place. As no H I is found in the 
inner regions, this atomic gas must have been converted into other 
forms. Indeed, hydrodynamical simulations have shown that in the 
rapidly changing potential of major mergers, the gas is unable to 
avoid collisions, thus dissipating energy and moving inward (e.g., 
Mihos & Hernquist 1994). The high-density clouds within this in- 
falling gas may well be compressed into molecular form (and sub- 
sequently form stars and/or star clusters, depending on the gas com- 
pression achieved (Jog & Solomon 1992; Elmegreen & Efremov 
1997)), while the low-density gas may shock-heat to the virial tem- 
perature (radiating in the soft X-rays). Indeed, NGC 7252 contains 
a large amount of cold molecular gas in its inner regions (Dupraz 
et al. 1990) as well as extended X-ray-emitting gas (Hibbard et al. 
1994). Radio continuum data of the central regions of NGC 7252 
(where the molecular gas is found, cf. Hibbard et al. 1994) indicates 
a modest star formation rate of 1 - 2 Mq yr -1 (Condon 1992). At 
this rate, the remaining molecular gas would turn into stars within 
2-4 Gyr (Dupraz et al. 1990). During that time, one would ex- 
pect the X-ray luminosity to grow as the remaining atomic and 
molecular gas contents are shock-heated or evaporated by elec- 
trons within the already existing hot plasma (e.g., Goudfrooij & 
Trinchieri 1998). Furthermore, if the merger produced a remnant 
with a deep enough potential well, the hot plasma would be sup- 
plemented by gas donated by mass loss from evolving giant stars 
and SNela and virialized in the potential well of the remnant (e.g., 
Brighenti & Mathews 1999). 

The observed properties of NGC 1316 fit in remarkably well with 
this evolutionary merger scenario: CO observations by Sage & Gal- 
letta (1993) revealed the presence of ~ 1 x 10 9 Mq of molecular 
gas. While this is a large amount relative to the molecular gas con- 
tent of "normal" giant elliptical galaxies (see, e.g., Knapp & Ru- 
pen 1996), the molecular gas is found to be already depleted in the 
central region, being concentrated in the dusty regions at ~ 30-50 
arcsec from the centre instead (Sage & Galletta 1993). In the con- 
text of the scenario presented above, this is entirely consistent with 
the situation expected a few Gyr after the merger took place. Fur- 
thermore, a recent analysis of optical spectra of the central 3" of 
NGC 1316 reveals a stellar population with a luminosity- weighted 
age of ~ 2 Gyr and a metallicity that is as high as other giant el- 
liptical galaxies of similar luminosity (Kuntschner 2000). This is 
again entirely consistent with the scenario of an evolved major 
merger: The merger occurred ~3-3.5 Gyr ago (as indicated by 
the age of the bright globular clusters) after which stars continued 
forming within the molecular gas in the central regions for a few 
Gyr (the latter process has now apparently largely ceased). In fact, 
the unusually small core radius and effective radius of NGC 1316 
(Schweizer 1980; Caon et al. 1994) are most probably due to this 
star formation episode in the central regions. Finally, the predicted 
build-up of hot, X-ray-emitting gas since the merger event has in- 
deed occurred in NGC 1316: analysis of ROSAT observations have 
shown the presence of ~ 10 9 Mq of hot (10 7 K) gas, part of which 
is (still?) distributed in tails and loops around the galaxy (Kim, Fab- 
biano & Mackie 1998). 

(ii) Tidal features are long-lived. Another important result of the 



numerical simulations of NGC 7252 by Hibbard & Mihos (1995) 
was that the tidal tails are not freely expanding. Instead, the most 
tightly bound bases of the tails were found to turn around most 
quickly, followed gradually by the material at larger distances from 
the centre. The rate of return was found to fall off roughly as t _5//3 , 
implying that tidal features such as those in NGC 7252 stay visi- 
ble for several Gyr. In the particular case of NGC 7252 they found 
that at a time of 4 Gyr after the merger^, roughly 20% of the cur- 
rent B-band luminosity of the tidal tails remains outside 4 effec- 
tive radii from the galaxy centre. It is therefore quite plausible that 
at least some of the large tidal tails and loops around NGC 1316 
still stem from a major merger event similar to the one that formed 
NGC 7252, i.e., one during which the 3 Gyr old, ~ solar-metallicity 
globular clusters were formed. However, we do not mean to imply 
that all tails and loops in NGC 1316 necessarily stem from one 
merger event. In particular, the tidal tail located ~ 6f 5 southwest of 
the centre (just outside our CCD field of view, but shown in Fig- 
ures 2 and 6 in Schweizer 1980, denoted Li) likely stems from a 
significantly more recent interaction with a small, gas-rich galaxy 
(cf. Mackie & Fabbiano 1998). 



6.3 Relevance to the nature of bimodal colour distributions 
in giant ellipticals 

Assuming that the luminous, metal-rich globular clusters in 
NGC 1316 were formed 3 Gyr ago during a major merger involving 
two spiral galaxies, what does the expected colour distribution of 
the resulting globular cluster system look like? Is it consistent with 
the observed one? Will it evolve into a bimodal colour distribution 
as seen in many "normal" giant ellipticals? 

and Figure ^| that the observed 



We have seen in Section 5.1 



colour distribution of genuine clusters in NGC 1316 is not bimodal 
but spans a rather narrow range around B — I= 1.85 (cf. also Grill- 
mair et al. 1999), consistent with the colour predicted by the BC96 
models for a 3 Gyr old, single-burst population of solar metallicity. 
In contrast, old, metal-poor, Milky Way halo globular clusters have 
B — I colours around 1.5-1.6. However, the post-merger cluster 
system of NGC 1316 should still contain many of the bulge/halo 
globular clusters that were associated with the progenitor galaxies. 
Is the observed colour-magnitude diagram consistent with the pres- 
ence of such "old" globulars as well? To address this question, we 
plot the predicted time evolution of the B — / colour and the abso- 
lute magnitude Mv for single-burst populations older than 1 Gyr 
in Figure [l2| (again using the BC96 models). For 14 Gyr old clus- 
ters, the observed B — I colour interval is found to be consistent 
with a metallicity range of 0.03 Z/Zq ^ 0.5, while the mean 
colour B—I = 1.85 would indicate Z ~ 0.1 Zq. While this seems 
somewhat higher than the median metallicity of the sample of halo 
globular clusters in our Galaxy (Z = 0.03 Zq, cf. Harris 19960), 
it is also somewhat lower than the median metallicity of the glob- 
ular cluster system that is associated with the bulge of our Galaxy 
(Z = 0.25 Zq; Minniti 1995). This actually suggests that the old 
clusters in NGC 1316 may be a mixture of clusters from two pro- 
genitor spiral galaxies with different bulge-to-disk ratios. To assess 
the reality of this suggestion, we encourage the development of nu- 
merical simulations of interactions of galaxies having bulge and 
halo globular cluster systems. The result of such simulations can 



§ Assuming H p = 75 km s 1 Mpc 1 

^ Compiled at http://www.physics.mcmaster.cayGlobular.htm 
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Figure 12. Time evolution of B — I colour index (top panel) and abso- 
lute magnitude My per 5 X 10 6 Mq (bottom panel) of single-burst stel- 
lar population models (Bruzual & Chariot 1996) for ages older than 1 Gyr. 
Model curves are plotted for a Salpeter (1955) IMF and the following metal- 
licities: 0.02 solar (short-dashed lines), 0.2 solar (long-dashed lines), 1.0 
solar (solid lines), and 2.5 solar (dot-long-dashed lines). Ages of 3 Gyr 
and 14 Gyr are indicated. The region hatched by dotted lines in the top 
panel depicts the B—I colour interval populated by genuine star clusters in 
NGC 1316 (cf. Figure ^|). The bottom panel shows the predicted evolution 
of V-band magnitudes (absolute magnitudes on the left-hand side, apparent 
magnitudes at the distance of NGC 1316 on the right-hand side) for a glob- 
ular cluster having the mass of lo Cen. The horizontal d otted lines indicate 
V magnitudes for specific cases, as discussed in Section |6~3| . 



be expected to significantly further our understanding of the nature 
of globular cluster systems with bimodal colour distributions, and 
their associated power as witness of the processes occurring during 
galaxy formation. 

Finally, the lower panel of Figure [l2| depicts the time evolu- 
tion of the absolute magnitude My (as well as V at the distance of 
NGC 1316) for a globular cluster similar to u Cen (M = 5 x 10 6 
Mq, Meylan et al. 1995). Taking this to be a typical brightest glob- 
ular cluster of any spiral galaxy, one can see that the "old", metal- 
poor clusters are expected to be fainter than V ~ 21.0. After the 
second-generation, solar-metallicity clusters evolve to an age of 14 
Gyr, the BC96 models predict a cluster like u) Cen to be V ~ 21.7 
and B—I ~ 2.2, i.e., those clusters will have AV ~ 0.7 mag fainter 
and A(B — I) ~ 0.3 mag redder (for a given cluster mass) than the 
clusters associated with the progenitor galaxies. These numbers are 
remarkably consistent with those of typical globular cluster systems 
(associated with giant ellipticals) which exhibit bimodal colour dis- 
tributions (e.g., Zepf, Ashman & Geisler 1995; Forbes et al. 1998; 
Grillmair et al. 1999). 

We conclude that the observed colour distribution of clusters 
in NGC 1316 is entirely consistent with that expected from a mix- 
ture of ( i) intermediate-age clusters of age ~ 3 Gyr having near- 
solar metallicities and (ii) "old" (~ 14 Gyr), metal-poor clusters 



that were associated with progenitor galaxies. As such, the colour 
distribution of this cluster system is expected to evolve into a bi- 
modal one, having properties consistent with those of well-studied 
giant ellipticals. The details of the appearance of the final colour 
distribution of the system will depend on the ratio of the number of 
newly formed clusters to the number of "old" clusters, which is lit- 
tle constrained at this point. We do know, however, that the massive 
end of the mass function of the newly formed clusters is at least as 
well populated as in the Milky Way and M 31 (i.e., as systems with 
several hundred clusters). 



7 MAIN CONCLUSIONS 

We have obtained spectra for 37 candidate star clusters in the field 
of the giant early- type galaxy NGC 1316, an obvious merger rem- 
nant featuring extensive shells, tails, and dusty filaments. Our main 
conclusions are as follows. 

• The radial velocity measurements reveal that 24 targets are 
genuine star clusters, and 13 targets are Galactic foreground stars. 
Almost all (22/24) genuine star clusters fall in the narrow colour 
range 1.6 & B — I 2.1, i.e., only 2 red targets with B — I > 2.1 
are highly reddened star clusters, the rest are foreground stars. 

• For the star cluster sample, we measure a mean heliocen- 
tric velocity Vhci — 1698 ± 46 km s~ and a velocity dispersion 
a = 227 ± 33 km s _1 within a galactocentric radius of 24 kpc. 
Partly responsible for the velocity dispersion is a significant rota- 
tion in the star cluster system, with a mean velocity of ~ 175 ± 70 
km s" 1 along a position angle of ~ 6° ± 18°. Using the projected 
mass estimator and assuming isotropic orbits, the estimated total 
mass is (6.6 ± 1.7) x 10 11 M within a radius of 24 kpc. The 
mass is uncertain by about a factor of two, depending on the or- 
bital assumptions and adopted distance. The above mass implies 
a Ml Lb ratio in the range 3-6, which is low for giant ellipti- 
cal galaxies. This may indicate that the stellar population of the 
integrated light of NGC 1316 is, in a luminosity-weighted sense, 
younger than usual. 

• For three bright star clusters in our sample, the signal-to- 
noise ratio is good enough to measure equivalent widths of Ha and 
the Ca II triplet absorption lines with confidence. Comparing those 
line strengths as well as measured broad-band colours with recent 
single-burst stellar population models (BC96), the metallicities of 
those 3 star clusters are found to be solar to within ± 0. 15 dex, and 
their most likely age is 3 Gyr. 

• We discuss the properties of the main body of NGC 1316 itself 
and conclude they are consistent with having hosted a major merger 
3 Gyr ago as well. The presence of intermediate-age globular clus- 
ters in NGC 1316 shows once again that globular clusters with near- 
solar metallicity do form during galactic mergers, and, moreover, 
that they can survive disruption processes taking place during the 
merger (e.g., dynamical friction, tidal disruption), as well as evap- 
oration. In this respect, NGC 1316 provides a hitherto "missing" 
evolutionary link between the recently studied young merger rem- 
nants of ~ 0.5 Gyr age such as NGC 3921 and NGC 7252 on one 
side, and "old" giant ellipticals featuring bimodal colour distribu- 
tions on the other side. 
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